Efficient and coherent excitation transfer across disordered molecular networks 
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We show that finite-size, disordered molecular networks can mediate highly efficient, coherent 
excitation transfer which is robust against ambient dephasing and associated with strong multi-site 
entanglement. Such optimal, random molecular conformations may explain efficient energy transfer 
in the photosynthetic FMO complex. 
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Recently, a vivid debate arose on the physical mecha- 
nisms underlying efficient transport in organic molecules. 
In particular, the excitation transfer from the photore- 
ceptor to the chemical reaction center in photosynthetic 
light harvesting complexes succeeds with astonishingly 
high transfer efficiency. Since many of these biological 
systems exhibit disorder and are coupled to noisy envi- 
ronments, arguably all models 1, 2] so far build on the 
fundamental hypothesis that disorder induces destructive 
interference in the coherent quantum evolution. The lat- 
ter, in turn, hinders transport [3, 4], what can only be 
overcome by added noise, such as to restore the classically 
diffusive behavior. However, this hypothesis is valid only 
in the thermodynamic limits i.e. for very large molecu- 
lar structures, while coherent transport across finite-size 
disordered samples is characterized by large fluctuations 
under configurational variations The photosyn- 

thetic FMO complex, as one of the most carefully studied 
examples fsj, is clearly a very finite molecular complex, 
far from the thermodynamic limit. New experimental 
data additionally provide clear evidence that excitation 
transfer is predominantly coherent even at room temper- 
ature [9l4ll|. on transient time scales (of the order of a 
few lOOfs) much shorter than the typical environment- 
induced decoherence times [l2l| . Hence, the cause of the 
observed transport efficiency must be rooted in general 
properties of coherent quantum dynamics on finite molec- 
ular networks, on time scales s/iorter than those on which 
environmental decoherence fully develops its detrimental 
influence |5-7]. As we will show, rare incidences of con- 
structive rather than destructive interference of transition 
amplitudes from the photoreceptor to the reaction center 
indeed do provide a possible explanation for these obser- 
vations, and enable strictly coherent transport efficiencies 
up to 100%. 

With some physical abstraction, a light harvesting 
complex can be viewed as a fully connected, disordered 
finite graph with N vertices. Under the assumption of 
purely coherent transport, it is the spatial distribution 
of these which fully controls the relative phases of the 
transition amplitudes which need to be added coherently 



to infer the transport efficiency from the input to the 
output site. Optimal transport efficiency is then equiva- 
lent to molecular conformations with strictly constructive 
interference of all these amplitudes - just conversely as 
in the case of disorder-induced destructive interference 
which dominates in the thermodynamic limit. On a fi- 
nite molecular network, this represents an optimization 
problem which can be solved by evolution, and suggests 
a statistical analysis, as follows. 

Coherent transport of a single excitation across a sam- 
ple of molecular sites is generated by the Hamiltonian 



H 



N 
i^3 = i 



(1) 



where and tri*'' mediate excitations and de- 

excitations of sites j and i from the local electronic 
ground state to the local excited state and vice versa, 
respectively. The excitation transfer aj^ a_ from site i 
to site j has a strength — Vj^i which depends on the 
specific nature of the inter-site coupling - that we assume 
to be of resonant (isotropic) dipole type, Vij = a/rfj, 



with 



and fj the position vectors of in- 



dividual sites. Input and output site define the poles 
of a sphere of diameter d which, via the coupling con- 
stant Win, out = ot/d'^, sets the natural time-scale of the 
dynamics induced by H . The positions of the remaining 
molecular sites are randomly (uniformly) chosen within 
this sphere, what induces a random distribution of the 
remaining Vi^j . 

To assess the probability for complete and rapid trans- 
fer of an excitation from the input to the output site we 
sample over different spatial configurations. Our figure 
of merit is the maximum probability - henceforth "trans- 
port efficiency" - 



pZI = max |(out|7/;(i))|' 
te[o,71 



(2) 



that an excitation injected at input is received at output 
after times no longer than T = 0.1 x 7r/(2|vin_out|)- This 
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FIG. 1. (Color online) a Fluctuation of the transfer efficiency pj^^t from input to output, for 500 different random conformations 
N = 7 sites (dots). Horizontal grey lines indicate the transfer efhciency of the experimentally infered [13] FMO Hamiltonian, 
as well as that of the optimal configuration FMO* compatible with the experimental error margin, b Probability densities 
-P(Pout ) of the transfer efficiency p^^t for 2.5 x 10* different conformations. For fully coherent dynamics (black curve) the 
mean value of p^^t amounts to 4.9%, and only 4.5 out of 10® configurations provide efficiencies larger than 90%. Under local 
dephasing (red and green curve) the mean efficiency drops to 3.9%. c Gains (green) and losses (red) of the transfer efficiency 
with dephasing. 



is one tenth of the time span the full excitation trans- 
fer would require if no intermediate sites were present. 
Note that this specific choice of T is immaterial for our 
subsequent discussion, provided T is sufficiently small as 
compared to the time scale set by Win,out j and long enough 

(T) 

to allow maximum values of Pou/ close to unity. This also 
implies (and is confirmed by inspection of the underly- 
ing numerical data) that large values of the transport 
efficiency as defined in © imply large values of the time- 
average of \{ont\ip{t))\^ over [0,7^ (an efficiency quanti- 
fier used elsewhere and vice versa. Only on time 
scales which are long with respect to T could both quan- 
tifiers lead to different predictions. On such time scales, 
however, coherent effects will fade out. 

Fig. [T^ shows the variation of the transport efficiency 
Pq^i- for a sample of 500 different random distributions of 
N = 7 sites as in the FMO complex fl^. p^J2 fluctuates 
wildly for different random configurations, as anticipated 




FIG. 2. (Color online) a Optimal spatial configuration of 
N = 7 sites offering fast, robust, and complete transport from 
input to output, b Time evolution of the on-site probabilites 
|{i|i/'(t))P generated by the Hamiltonian defined by a. i is ei- 
ther the input site , the output node , or an interme- 
diate site . At time T', only the output site is populated. 

At intermediate times t < T' , the excitation is spread over 
several sites, leading to high values of the bipartite and 
quadripartite entanglement, see text. 



in our motivation above. Remarkably, very high trans- 
port efficiencies can be achieved as evident from Fig. [TJd, 

(T) 

where the probability density of p^^l obtained for a sam- 
ple of 2.5 X 10® realizations is shown in black. Efficien- 
cies above 90% are achieved for approx. five configura- 
tions out of a million, despite the fact that the average 
efficiency is only 4.5%. Therefore, evolution can choose 
from truly exceptional, albeit statistically unlikely molec- 
ular conformations with excellent transport properties. 
Such a configuration is depicted in Fig. [2] where - in con- 
trast to the random choice of conformations in Fig. [1] - 
we have maximized Pou{ by iterative optimization of the 
sites' positions. Thereby, we find certain conformations 
which achieve even 100% transport efficiency. The exam- 
ple shown in Fig. [2^ spells out that optimal arrangements 
are asymmetric and non-periodic, hence far from trivial 
{e.g., lattice-like, collinear) structures. 

To assess the efficiency of the actual FMO complex, 
in comparison to our present results, we employ an ap- 
proximate Hamiltonian [13\ inferred from experimental 
data, and obtain a transport efficiency of only 5.7% in a 
time window of 1.6 x 10~^'^s, hence close to the average 
value of our random model. The time window was here 
defined by the coupling strength between chromophores 
1 and 3, T — 0.1 x 7r/(2|wi 3|). However, variation of 
the off-diagonal (diagonal) matrix elements by at most 
3.2 X 10"te"i (17 X 10"/is"i) - what is the absolute 
error margin deduced from the experimental data [l3| 
- is compatible with an alternative, optimal configura- 
tion (FMO*) with a transport efficiency of 43.1%, see 
Fig. [Ih {h is Planck's constant)! Furthermore, much as 
the optimal configuration depicted in Fig. [2^, this config- 
uration's efficiency is robust under statistical variations 
with a spread of lO^te^^ (5.4 x lO^/is"!) on the FMO* 
Hamiltonian's off-diagonal (diagonal) elements, in the 
sense that such variation yields a Gaussian distribution 
43.1 ±5% of the transport efficiency around the optimum. 

Let us now consider the same transport problem in the 
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presence of environmental noise. Fig. [T): shows the effi- 
ciencies Pq^j for the same statistical sample as in Fig.[T]a, 

under local dephasing with a strong rate 7 = 2/T. Cases 
(T) 

in which p]^^l is decreased by dephasing are highlighted 
in red, enhancements of transport efficiency is depicted 
in green. The plot very neatly spells out a clearly di- 
chotomous impact of the environment: 

(i) Whenever constructive quantum interference en- 
hances transport in the absence of environment coupling, 

(T) 

the noise reduces the transport efficiency p^^l very con- 
siderably. 

(ii) In contrast, if quantum coherence suppresses trans- 
port in the strictly coherent case, dephasing will enhance 

(T) 

Pqjjj , though only very marginally so. 

(iii) Notwithstanding, even in the presence of the 
rather strong dephasing chosen for our simulation, those 
rare molecular conformations which provide efficient ex- 
citation transfer maintain this characteristic property un- 
der environmental coupling, just at reduced levels, and 
are still clearly distinct from those conformations which 
hinder transport. 

The crossing of the probability densities obtained with 
and without noise, see Fig.[Tj3, identifies a level of 7.6% as 
the demarcation line between transport efSciencies which 
are strongly reduced (Pout > 7.6%) or marginally en- 

(T) 

hanced {p'J^l < 7.6%) by added noise. This is again in 
qualitative accord with the reported data on the FMO 
Hamiltonian, as well as with its efficient variant FMO* 
introduced above: while the FMO efficiency increases 
from 5.7 to 12.3 percent, the FMO* efficiency is reduced 
from 43.1 to 21.5%, in the presence of noise. 

Note that such dichotomous behavior as identified here 
for finite systems on transient time scales is absent in 
the thermodynamic limit of infinite systems and/or suf- 
ficiently long transport times, where noise completely 
destroys quantum coherences and tends to induce near- 
classical, diffusion-like behavior. It is well-known that 
noise then enhances transport by suppressing destructive 
quantum interference [l|, , though can in general not 
compete with the transport efficiencies brought about by 
constructive quantum interference - even when the envi- 
ronment coupling strengths are optimized 

Efficient quantum transport as observed above relies 
on the constructive interference of a large number of (in- 
put to output) transition amplitudes, reminiscent of ef- 
ficient quantum algorithms. Therefore, let us now ad- 
dress the question whether multi-site entanglement is of 
similar relevance for the molecular transport problem as 
it is for quantum computation - an issue of much re- 
cent interest 14, [l^. In our present problem, pre- 
cisely one excitation propagates from input to output, 
so that the transporting states are close relatives of the 
W-states [3] - a well-known class of entangled states 
of multipartite systems, that are known to be particu- 
larly robust against decoherence induced by dephasing 



or spontaneous decay 17|. These states also provide a 
clear relation between the excitation's localization and 
the system's entanglement properties: The latter can 
be characterized in terms of a hierarchy of quantities 
{c,j{ip)}^=2,...,N which are strictly positive if entangle- 
ment in \tp) is shared by at least v sites, and vanish oth- 
erwise. In particular, each of the c^{ip) is a function of 
the statistical moments Mk{ip) = J2jLi IOIV')P*'': where 
{\j)}j=i,....N is the canonical site basis. In terms of the 
second moment M2{tp) - which is nothing but the in- 
verse participation ratio f3| frequently used in statistical 
descriptions of complex quantum systems - C2 (■(/') (which 
is a multi-partite generalization [ISj of a standard bipar- 
tite entanglement measure [l^) reads 



C2 i^) = V^T^(1-M2("^ . (3a) 

Higher order c^{'ip) are analogous functions of the mo- 
ments Mi{tp) = 1 to M^{ip). For instance, the quadri- 
partite measure reads 



C4(V') 



1 - 6A/2 + 8M3 + 3Af| - 6M4 \ ^'"^ 
1 -6/7V+ ll/iV2 -6/iV3 j 



(3b) 



Under purely coherent dynamics, multi-partite entangle- 
ment thus encodes the detailed localization properties of 
the excitonic wave function. 

With these tools at hand, we can now correlate the 
transport efficiency with the multi-site entanglement 
which is generated during the transport process. Figs.|3h 
and [5)d show the probability density of the transport ef- 



ficiency, P(piut^ 



sr') 



and P(p^J2 



conditioned 

(T') 



on the maximal bi- and quadripartite entanglement, C2 



and 



Sr') 



which is generated during the exciton propa- 



gation from input to output: 



c^J ^ = max Cy{i>{t)) 
te[o,r'] 



(4) 



T' < 7" is the time at which the maximum output prob- 
ability p^J^ is reached, see Fig. [2Jd. Clearly, efficient 
transport necessarily requires strong entanglement. This 
is most prominently spelled out in Fig. [3^, where high 
transport efficiency {e.g. pl^^l > 0.5) is only reached at 

(T' ) 

high values of the entanglement (cj > 0.8) shared be- 
tween at least two sites. Note the kink in the distribution 
visible at = •\/7/12 ~ 0.76, which corresponds to 

maximal entanglement between exactly two of iV = 7 
sites. Here, transport is inhibited since the excitation 
may be trapped in a singlet state between two sites which 
are accidentally placed very close to each other. 

The correlation between transport and entanglement 
visible in Fig. [3^ prevails for higher orders of the c'iT \ 
but is less pronounced for increasing v, as evident from 
the exemplary case of c^J^ ^ in Fig. IJJj. Here, although 
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FIG. 3. a Contour plot of the conditional probability density P{p'^^ j i.e. the probability distribution of the transport 

efficiency p'^^l across 7 sites, given the generation of a certain maximal level Cj of at least bipartite entanglement, during the 

transfer time T'. b Same as a, but for quadripartite entanglement c^^ \ c Same conditional probability distribution as in a, 
when all sites are locally coupled to a dephasing environment, with dephasing rate 7 = 2/T. In all three cases, large transport 
efficiencies require a minimum amount of entanglement. 



(T) 



0.2 are pos- 



moderate transport efficiencies like p, 
sible at very small values of , fiigher transport ef- 
ficiency still requires a certain amount of entanglement 
(i.e. C4 > 0.5) between at least four sites. 

Finally, to gauge the robustness of the observed corre- 
lation under decoherence, Fig. [3t shovifs the conditional 



probability density P(p, 



(T) I AT')\ 



estimated for mixed 



states \2Ql when the individual molecular sites are locally 
coupled to dephasing environments, with the same decay 
rate 7 = 2/T as in Fig.[T]:. Fully consistent with our dis- 
cussion of Fig. [1] the correlation between entanglement 
and transport efficiency remains qualitatively unaffected, 
however with smaller transport efficiencies and entangle- 
ment levels than in the strictly coherent case. This cor- 
relation is particularly remarkable here, since the above 
equivalence of multi-site coherence and multi-site entan- 
glement cannot be established any more under open sys- 
tem dynamics. 

In conclusion, we have seen that very well defined 
molecular configurations, which can be found by iter- 
ative optimization, mediate highly efficient and robust 
transport across molecular networks alike the FMO en- 
ergy harvesting complex. Even in the presence of rather 
strong dephasing does efficient excitation transfer due to 
constructive quantum interference remain a distinctive 
feature of these conformations. Efficient transport is fur- 
thermore conditioned on the build-up of strong inter-site 
entanglement in the course of the exciton transfer. This 
is clear evidence of the functional role of multi-site entan- 
glement on the level of biomolecular (quantum) dynam- 
ics. Whether, beyond that, biology has ways to harvest 
the statistical, non-local quantum correlations between 
single excitation events at different sites of VF-like states 
on FMO-like functional units remains an intriguing ques- 
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